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Abstract 

We examine two semi-analytical methods for estimating the baryon asym- 
metry of the universe (BAU) generated in scenarios of "local" electroweak 
baryogenesis (in which the requisite baryon number violation and CV violation 
occur together in space and time). We work with the standard electroweak 
theory augmented by the addition of a CV violating dimension six operator. 
We work in the context of a first order phase transition, but the processes 
we describe can also occur during the evolution of a network of topological 
defects. Both the approaches we explore deal with circumstances where the 
bubble walls which convert the high temperature phase to the low tempera- 
ture phase are thin and rapidly moving. We first consider the dynamics of 
localized configurations with winding number one which remain in the broken 
phase immediately after the bubble wall has passed. Their subsequent decay 
can anomalously produce fermions. In a prelude to our analysis of this effect, 
we demonstrate how to define the C and CV symmetries in the bosonic sector 
of the electroweak theory when configurations with nonzero winding are taken 
into account. Second, we consider the effect of the passage of the wall itself on 
configurations which happen to be near the crest of the ridge between vacua 
as the wall arrives. We find that neither of the simple approaches followed 
here can be pushed far enough to obtain a convincing estimate of the BAU 
which is produced. A large scale numerical treatment seems necessary. 
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I. INTRODUCTION 



One of the central challenges of modern cosmology is to explain how the baryon asym- 
metry of the universe (BAU) can be dynamically generated from baryon number sym- 
metric initial conditions. The relevant observational quantity is the baryon to entropy 
ratio which is tightly constrained by primordial nucleosynthesis to lie in the interval 
2 x 10~ 10 < ub/s < 5 x 10~ 10 . Sakharov JT] showed that there are three necessary condi- 
tions which must be satisfied if a particle physics model is to produce a net cosmological 
baryon asymmetry. These are a violation of baryon number (£>) conservation, a departure 
from thermal equilibrium, and violations of charge conjugation (C) and charge-parity (CV) 
symmetries. The standard model of electroweak interactions includes all three ingredients, 
as we discuss below. This observation has led to the suggestion that the standard model 
might be responsible for the generation of the baryon asymmetry of the universe @-16|. 

Baryon number is exactly conserved at the classical level in the standard electroweak 
theory, but at the quantum level this ceases to be true as a consequence of the anomaly [fTT 



The relevant process at zero temperature is an exponentially suppressed tunnelling under 
the energy barrier which separates inequivalent vacua in gauge and Higgs field configuration 
space. The lowest point on the crest of the barrier is a saddle point configuration called the 
sphaleron [|T^]. At nonzero temperature and in particular for temperatures around or above 
the critical temperature of the electroweak phase transition it is possible for thermal effects 
to cause classical transitions over the sphaleron barrier [pl|l9|- Although at temperatures 
above the electroweak phase transition anomalous baryon number and lepton number (£) 
violation is unsuppressed in the standard model, this occurs in such a way that B — £ 
is conserved. Thus, any preexisting B — £ asymmetry cannot be erased by electroweak 
processes. We are assuming that above the temperature of the electroweak phase transition 
B and £ are both zero. 

In order to construct a scenario for electroweak baryogenesis we must achieve a departure 
from thermal equilibrium. In the context of the early universe this is typically satisfied in one 
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of two ways. One possibility is that the electroweak phase transition is considered to be first 
order so that the violent conversion of the high temperature phase to the low temperature 
phase results in non-equilibrium conditions near the bubble walls separating the phases 
f5HT4| . An alternative is that non-equilibrium processes occur in topological defects, such as 
cosmic strings, around which the electroweak symmetry is restored [|15|| . In defect scenarios, 
the phase boundaries between the interior and exterior of the defects behave in an analogous 
manner to bubble walls. In this paper, we shall phrase our considerations entirely in terms 
of bubble walls at a first order phase transition, even though the processes we discuss also 
arise in defect mediated baryogenesis. 

We make several standard assumptions related to the rate per unit volume T for baryon 
number violating processes in thermal equilibrium, which we write as 

T = K(a w T) 4 . (1.1) 

We assume that the phase transition is strongly enough first order that at temperatures below 
T c , T is small enough that the BAU created at the transition is not erased by subsequent 
electroweak baryon number violating processes. (The constraint turns out to be k < 10 -9 
for T < T c PC| , PT]] .) It now seems likely that this requires extending the field content of 
the theory in some way which makes the phase transition more strongly first order. As is 
often done, we use the minimal standard model with a Higgs mass lighter than is allowed 
by experiment to model the strongly first order phase transition in the extended theory. 
We also assume, again as is standard, that at temperatures just above T c the dimensionless 
quantity k is not very small. 

Sakharov's third condition is that C and CP must be violated. The standard electroweak 
theory is maximally C violating due to the V-A nature of the interactions. As we will show, 
the purely bosonic sector of the theory is also C violating once topological properties of 
the sector are treated correctly. The only violation of CP in the electroweak sector of the 
standard model occurs in the CKM matrix. Electroweak baryogenesis scenarios using this 
source of CV violation are generally expected to produce a BAU which is far too small. We 



can also ignore the CP violating 9 term of QCD since this is known to be tiny. Therefore, if 
we are to generate an appreciable baryon asymmetry it is necessary to consider extending 
the standard model to include new sources of CP violation. A popular approach has been 
to consider a two-Higgs doublet model with explicit renormalizable CP violating terms in 
the Higgs sector P,p|,P, |12| , p!3|,[15| . This enhances the CP violation beyond that present in 



the standard model but adds a host of new free parameters to the theory. Here we follow 
the simpler approach first discussed by Dine et al. [0. We assume that there is some CP 
violating physics (beyond that in the standard model) at an energy scale M greater than 
v = 250 GeV. At scales less than M, the effective theory is the standard model plus 
nonrenormalizable operators, some of which are CP odd. In this paper, for definiteness, we 
add the dimension six CP odd operator 

= ^Tr($t$)Tr(F,^) , (1.2) 

to the standard model Lagrangian density, and do not augment the field content of the 
theory. Here $ is the Higgs field defined in (j2.3|), F^ u is the SU{2) field strength tensor 
defined in ( |2.2|) and 

pr> = : l e ^F a p . (1.3) 

is the dual of the field strength tensor. O is the lowest dimension CP odd operator which 
can be constructed from minimal standard model Higgs and gauge fields. Standard model 
CP violation in the CKM matrix does induce the term O with M = v in the effective action, 
but the coefficient b is thought to be tiny. For us, b is a dimensionless coupling measuring 
the strength of the CP violating physics at the scale M, above which our effective theory 
ceases to be valid. Throughout most of this paper, b and M occur only in the combination 
b/M 2 , and in this sense we are introducing a single new parameter. 

The operator O induces electric dipole moments for the electron and the neutron, and 
the strongest experimental constraint on the size of such an operator come from the fact 
that such dipole moments have not been observed. Working to lowest order (one-loop) we 
find 



d e _m e sm (9 W ) b ^ ( M 2 + m 2 H 

f2 



, y2 lH— T^l- (1-4) 



M arises in the logarithm without b because M, the scale above which the effective theory 
is not valid, is the ultraviolet cutoff for the divergent loop integral. A result similar to ( |1.4|) 
was obtained in Ref. Using the experimental limit d e /e < 4 ■ 10 _27 cm we find the 
bound 
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The experimental limit |24| on the neutron electric dipole moment d n is weaker than that on 
d e , but because d n is proportional to the quark mass rather than to the electron mass, the 
constraint obtained using d n is comparable to (|1.5| ). A baryogenesis scenario which relies on 
CV violation introduced via the operator O must respect the bound ( |1.5| ). 

We wish to estimate the BAU generated during a first order electroweak phase transition 
in the minimal standard model with the addition of the operator O with a coefficient b/M 2 
satisfying ( |1.5|) . Baryogenesis can occur either "locally" (baryons are produced as a result 
of B violating processes and CV violating processes occuring together near bubble walls) or 
"nonlocally" (baryons are produced as a result of CV violation occurring as particles bounce 
off bubble walls and B violation occuring away from bubble walls in the high temperature 
phase). In general, both local and nonlocal baryogenesis will occur and the BAU will be 
the sum of that generated by the two processes. We have nothing to add to the existing 
treatments P,|9Hl4| of nonlocal baryogenesis, also called baryogenesis by charge transport. In 



this paper, we reconsider the models of baryogenesis first introduced by Turok and Zadrozny 
H and Dine et al. , and try to estimate the BAU produced locally. Throughout this paper, 
we assume that the time during which the expectation value of the Higgs field is changing 
at any point in space is short compared to all other timescales in the problem. Throughout, 
we refer to walls satisfying this criterion as "thin" , even though they may meet the criterion 
either by being thin enough or by moving fast enough. It is worth noting that for walls which 
move at a supersonic velocity, diffusion from the wall foward into the high temperature phase, 



and hence nonlocal baryogenesis, are not possible. In the opposite limit to that we treat, 
that is if the walls are "thick" in the sense that the Higgs field changes slowly relative to 
all other timescales, then the effects of the operator O can be treated approximately as a 
chemical potential for baryon number 0, realizing a possibility first considered in Ref. Q 
and called spontaneous baryogenesis. The local contribution to the BAU in the thick wall 
limit was first estimated by Dine et al. [0] and has been further analyzed in Refs. [p ,1X1,121, 
and is not affected by the considerations of this paper. 

In this paper we examine two semi-analytical methods for estimating the efficiency of 
local baryogenesis occurring in the standard model with the addition of the operator O, 
assuming a strongly first order phase transition and thin bubble walls. In section III, we 
explore an approach pioneered by Turok and Zadrozny || to estimate the baryon asymmetry 
by considering the relaxation of topologically nontrivial field configurations produced during 
the phase transition. In contrast with Ref. || we study the full Higgs plus gauge dynamics 
and, because we introduce CV violation via the operator O rather than using the two Higgs 
doublet model employed in Ref. [[J, we need not and do not add new fields to the theory. 
Turok and Zadrozny considered effects which occur while a (thick) bubble wall is passing. 
We use their method to consider physics after the passage of a thin wall. Unfortunately, we 
show that the use of a one-parameter family of spherically symmetric configurations as in 
Ref. || can give qualitatively misleading results, and demonstrate the difficulty of obtaining 
analytic estimates in this picture without doing a full scale numerical simulation. In Section 
II we discuss certain properties of gauge and Higgs field configurations in the standard 
electroweak theory, as a prelude to our discussion of baryon number production in Section 
III. Included in Section II is a demonstration that when the C and V transformations are 
defined to treat configurations with winding properly, the bosonic sector of the electroweak 
lagrangian is not C or V invariant but is CP invariant. In most treatments, C violation 
arises because of the V — A nature of the fermionic part of the theory, so it is amusing to 
see that even in the bosonic sector, C is not a good symmetry. In Section IV we turn to 
a method introduced by Dine et al. J7j. We consider configurations which happen to be 



near the crest of the ridge between vacua as the wall arrives, and estimate the extent to 
which their velocity in configuration space is modified by the operator O during the passage 
of the wall. If the wall is thick, it turns out that the velocity of motion in configuration 
space is not affected. The asymmetry arises because O affects the potential energy surface 
in configuration space during the passage of the wall. If the wall is thin, as we assume in 
this paper, there is no significant time during which the potential energy surface is affected, 
but the configuration space velocities are affected asymmetrically. We show how to estimate 
the quantities entering the final result in the thin wall limit correctly. Nevertheless, we 
argue that difficulties of the kind encountered in Section III also apply to the method of 
Section IV, rendering the estimate for the BAU more of an upper bound than an estimate. 
Thus, it seems to us that neither approach can be pushed far enough to obtain a convincing 
semi-analytical estimate of the BAU produced by local electroweak baryogenesis in the thin 
wall limit, and a large scale numerical simulation is called for. 

II. RELEVANT PROPERTIES OF THE STANDARD ELECTROWEAK THEORY 

In this section we examine certain properties of the standard model which are important 
for the arguments of Section III. Although our motivation is to set the stage for the next 
section, we have also endeavoured to write this section in such a way that it is independent of 
the rest of this paper. The relevant dynamics take place in the bosonic sector. We consider 
fermion production in the background of the evolving Higgs and gauge fields. We do not 
take into account the back-reaction of the fermions on the bosonic fields. Thus, we begin 
by looking at the purely bosonic part of the standard model where, for simplicity, we ignore 
the U(l) hypercharge gauge field. 



C = -^Tr(F^F^) 




— v 




where 




dpA v - d u A^ - ig[A^ A v ] 




(<9 M - igA^Q 



(2.2) 
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with = A a ^r a /2 where r a are the three Pauli matrices. The standard Higgs doublet 
if = (</?!, ip 2 ) is related to the matrix <3> by 

*<*•«> = (-% %)■ (2 ' 3 » 

Here v = 247GeV and g = 0.65. The gauge boson mass is m w = \gv and the Higgs boson 
mass is ran = V2Xv. 
Note that 

&<S> = (<p* m + ip* 2 <p 2 ) (J J) , (2.4) 

so that we can write 

*=%U, (2.5) 

where a 2 = 2 (</?*</?i + ^2^2) = Tr$t<I>, and [/ is an SU(2) valued field which is uniquely 
defined at any spacetime point where a does not vanish. Without loss of generality we 
impose the condition that at all times 

lim <r(x,t) =v , (2.6) 

|x|— >oo 

,^^')=(J ?) • < 2 - 7 > 

In A = gauge, a vacuum configuration is of the form 

v 

$ = -= U 
V2 

A a = —djUU* . (2.8) 

At any time t when er(x, t) 7^ for all x we have that U (x, t) is a map from R 3 with the 
points at infinity identified, that is S* 3 , into SU(2) and therefore C/(x, t) can be associated 
with an integer-valued winding 

N H (t) = w[U] = / ^e^'TrfC/t^C/C/tg,.^^^^ ? (2 .9) 



the Higgs winding number. If $(x, t) evolves continuously in t then Nn(t) can change only 
at times when there is a zero of a at some point in space. At such times, Nh is not defined; 
at all other times, it is integer-valued. Note that the Higgs winding number of a vacuum 
configuration (|2.8| ) is equal to its Chern-Simons number 

N CS (t) = g|^/ d 3 xe^ k Tr U^Aj, + ^igA^A^ . (2.10) 

For a general non-vacuum configuration the Chern-Simons number is not integer-valued. 

A. Topologically Interesting Configurations 

In this section we are interested in the dynamics of nonzero energy configurations with 
nonzero Higgs winding. A simple example is 

$(x) = ^[i](x) 
A„(x) = , (2.11) 

where £/[i](x) is a winding number one map, say, 

C/m(x) = exp (^(r)r ■ x) , (2.12) 



with 77(0) = — 7r and 77(00) = 0. The configuration ( |2.11| ) has no potential energy but does 



carry gradient energy because the covariant derivatives D$> do not vanish. This configura- 



tion has Nh = 1. If the configuration (|2.11| ) were released from rest it would radiate away 
its energy and relax towards a vacuum configuration. There are two very different ways for 
this to occur 0. If the characteristic size of C/[i] is large compared to m^ 1 , then the gauge 
field will evolve until it lines up with the Higgs field making the covariant derivatives zero, 
and at late times Nh will still be one. If the characteristic size is small the configuration 
will shrink, the Higgs field a will go through a zero, and at late times Nh will be zero. This 
dynamics is the subject of the next section. 

Note that Nh is not invariant under large gauge transformations. However, the change 
in Higgs winding, ANh, is gauge invariant and the two distinct relaxation processes are 
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distinguished by whether ANh is zero or nonzero. Throughout this section we choose the 
gauge such that our prototypical initial configuration is of the form ( 2.11|) which has Nh = 1. 



The configuration Q2.11| ) is similar to the Skyrmion which is a winding number one 
soliton in the nonlinear sigma model associated with QCD. The Skyrme lagrangian has a 
four derivative interaction which is not present in ( |2.1|) . This term stabilizes the soliton 
because it prevents the winding number one configurations from shrinking to zero size. 
Note that in the Skyrme model winding number one solitons are identified with baryons and 
winding number minus one solitons are viewed as antibaryons. 

We could view the action ( [2.1|) as an effective theory which describes the low energy 



degrees of freedom of some more fundamental theory such as technicolor. The sigma model 
of ( |2.1| ) would be to technicolor what the usual sigma model which describes pions is to 
QCD. The addition of a sufficiently large Skyrme-like term would result in classically stable 
electroweak solitons which we could call techni-skyrmions. The Higgs winding of such an 
electroweak soliton would be identified with its technibaryon number. In this paper we are 
not modifying ( |2.1| ) except for the addition of the dimension six CP violating term fll.2| ) and 
all configurations with winding are unstable. 

B. Fermion Production 

A wound up configuration of the form Q2.11| ) is not stable and if released from rest it will 
evolve to a vacuum configuration of the form ( [2.8|) plus radiation. In the process fermions 
may be anomalously produced. If the fields relax to the vacuum by changing the Higgs 
winding then there is no anomalous fermion number production. However, if there is no net 
change in Higgs winding during the evolution (for example a never vanishes) then there is 
anomalous fermion number production. 

To understand these claims consider two sequences of configurations beginning with the 
wound up configuration (|2.11| ) and ending at the classical vacuum ( |2.8|) . The first sequence 



ends at the vacuum ( |2.8| ) with U = 1 while the second ends up at U = Uny Note that 



these sequences cannot be solutions to the classical equations of motion since the initial 



configurations carry energy whereas the final ones do not. Throughout both sequences we 
maintain the boundary conditions ( |2.6| ) and ( |2.7| ). For the first sequence, a must vanish at 
some intermediate configuration since the Higgs winding changes. For the second sequence, 
the change in Higgs winding is zero and a need not vanish. 

Now introduce an SU{2)i weak fermionic doublet, ip. The fermion is given mass through 
the usual gauge invariant coupling to the Higgs field $ and for simplicity we assume that 
both the up and down components of the doublet have the same mass, m. The fermion field 
is quantized in the background of the bosonic fields given by our interpolation. 

Now, the anomaly equation 

^ = J^Tr(FF) , (2.13) 

when integrated, implies that the change in the fermion number from the beginning to the 
end of a sequence is given by 

-w[U) , (2.14) 

where U is that of the final configuration Q2.8| ). For the first sequence w is one whereas for 
the second it is zero. Thus fermion number is violated in processes for which the configu- 
ration ( p.ll| ) unwinds via gauge unwinding, but is not violated when such a configuration 
unwinds via a Higgs unwinding. 

For both of the interpolations which we have considered, the final background configu- 
ration is a vacuum configuration of the form ( |2.8| ). In this background the lowest energy 
fermion state has fermion number zero and the fermion number of any other state with n ip 
particles and m if) antiparticles is n — m. However, the fermion number of the initial state, 
where the background has the non- vacuum form ( |2.11| ), is more complicated. Suppose that 



f d 3 xJ° 




fd 3 xJ° 






final 




initial 



the winding number one map Um in ( |2.11| ) has a characteristic size L. It is known that 



if mL ^> 1 then the lowest energy fermion state in the presence of this background has 
fermion number one [P5|-|2"T|]. In this sense the configuration Q2.11D is said to carry fermion 



number. If our initial fermion state is this lowest energy state and we let this state evolve 
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in the changing background of our first interpolation then the final state will contain one 
net ip particle. Although a ip particle must be produced, there is no violation of fermion 
number since both initial and final states have fermion number one. If we evolve the same 
initial state in the background of the second interpolation there will be no net ip particles in 
the final state. This is consistent with anomalous fermion production as described by equa- 
tion (|2.14j) . On the other hand, if we begin with the configuration fl2.11|) and mL <C 1, then 
the lowest energy state in this background has fermion number zero. If this state evolves 
in the background of the first interpolation then there are no net ip particles in the final 
state whereas if it evolves in the background of the second interpolation there is one net 
antifermion produced. This discussion of where the fermion number resides does not alter 
the general conclusion which is that if we smoothly interpolate from fl2.ll ) to ( |2.8| ) with 



U — 1 then there is no anomalous fermion production in this background whereas if we end 
up at ( |2.8|) with U — Um then there is anomalous fermion production. 

Of course, we are actually interested in the dynamical evolution of configurations such 



as ( 2.11|) which are released from rest and end up as outgoing radiation. In this case it is 



dangerous to use the anomaly equation ( |2.13| ) since / d 4 x Tt(FF) is not well-defined as an 



integral and any answer can be obtained for the change in fermion number [EH] . Nonetheless, 



the results of our previous discussion still apply. If the configuration (|2.11| ) is released and 
falls apart without ever going through a zero of the Higgs field, then the analysis of Ref. |2!| 
is directly applicable and we conclude that one net antifermion is produced just as we did 
with the second of our interpolations. If the Higgs field unwinds by going through a zero, 
then one can use arguments presented in Refs. |2^,[29|] to demonstrate that the presence of 



outgoing radiation in the final configuration does not affect the result above, namely that 
there is no fermion number violation. 

C. Discrete Symmetries: C, V and CV 



In this subsection we study the properties of the lagrangian ( |2.1| ) under the discrete 
transformations C, V and CV. We must keep in mind that certain configurations of $ 

11 



can be associated with a particle number which equals the Higgs winding defined by ( |2.9| ) 
and ( |2.5| ). With the identification of winding with particle number, we will see that ( |2.1| ) is 
not invariant under C and V separately but is invariant under CP . 

First let us describe the parity transformation. Let V be the obvious discrete transfor- 
mation defined as 

V : 4,(x,t) ->Ao(-x,t) 

T 3 ' : $(x,t) -> $(-x,t) . (2.15) 

Under the parity transformation V' the lagrangian ( |2.1|) is invariant. However, under V' the 
winding number given by fl2.9| ) changes sign and we certainly do not want particle number 
to change sign under parity! Now parametrize the unitary matrix U appearing in Q2.5| ) by 
the fields 7r a (x, t) which are weak scale sigma model analogs of the low energy pions 

$(x,t) = ^l e xp(W(x,t)) . (2.16) 

We can define a parity operator V which takes 7r a (x, t) to — 7r a (— x, t) which is the conven- 
tional transformation property of the ordinary pions. Consistent with this we define 

V : A)(x,t) -> A (-x,t) 

V : Aifat) -A(-x,t) 

P : $(x,t) -> $ f (-x,t) . (2.17) 

With this definition of parity, which we adopt, the winding w[U] is unchanged by a parity 
transformation. However, in this case the weak interaction lagrangian ( |2.1| ) is not parity 
invariant. 

Now turn to charge conjugation. If we include the U(l) interaction we see that A 3 ^ is a 
linear combination of the photon and the Z-boson. Since the photon is charge conjugation 
odd we certainly want — > — A z . Similarly we want W+ — ► — W~ and — > — 
which is the same as A 1 — > —A 1 ^ and A 2 ^ —>■ A 2 . This is equivalent to the requirement that 
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the matrix = A a ^r a /2 transforms into tiA^t^. Thus we can attempt to define charge 
conjugation by C where 

C : A^ -> r 2 A M r 2 

C':$^t 2 $t 2 (2.18) 

This transformation has the property that (fi —>■ ip{ and ^ — ► ^2 which are the expected 
transformations of the complex fields and furthermore leaves the lagrangian ( |2.1|) invariant. 
Note that the Lagrangian (|2.1 ) is invariant under the transformation $ — > B^B, A^ — > 



B^A^B, where 5 is a spacetime independent SU(2) matrix. With B = ir 2 we obtain the 



transformation (|2.18|) . Thus, the transformation C is one element of a continuous global 
symmetry group and should not be viewed as a discrete transformation. In addition, under 
C the winding number is unchanged which is unacceptable. 

Return to the parametrization of $ given by ( |2.16|) . Under charge conjugation we expect 
(7T 1 , 7T 2 , 7T 3 ) — > (7T 1 , — tt 2 , 7T 3 ). To be consistent with this let us define 

C : Af, -> r 2 A M r 2 

C:$^ T2 $V 2 (2.19) 

Adopting this definition we have that the winding flips sign under charge conjugation. How- 
ever, the weak interaction lagrangian (|2.1| ) is not C invariant. 
Consider CV defined as the composition of (|2.17|) and (|2.19|) : 



CV : A (x,t) — > T 2 Ao(-x,t)T 2 
CV : Ai(x,t) -> -r 2 A(-x,t)r 2 

CV : $(x, t) -> r 2 $(-x, t)r 2 . (2.20) 

All fields transform as expected, for example y?i(x, t) — > </?*(— x,t), the lagrangian ( |2.1| ) is 
invariant under CP and the Higgs winding changes sign. Note that CV = CV. In defining 
C and V one must be careful to ensure that winding flips under C but not under V. In so 
doing, one discovers that with appropriate definitions of C and V the bosonic sector of the 
electroweak lagrangian is not C or V invariant but is CV invariant. 
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III. LOCAL BARYOGENESIS THROUGH UNWINDING 



In this section we explore the possibility of using the method of Turok and Zadrozny || 
to estimate the baryon asymmetry produced by local baryogenesis in a scenario in which 
the electroweak phase transition is strongly first order and the bubble walls are thin. The 
theory we treat is the standard model augmented by the CP violating operator O of ( |1.2| ). 
Turok and Zadrozny II studied the classical dynamics of topologically nontrivial gauge and 
Higgs field configurations in the presence of CP violation. Following Ref. M, we begin by 
considering spherically symmetric nonvacuum configurations of the form ( |2.11| ) with Higgs 
winding Nh = ±1 and discuss their dynamics when they are released from rest and evolve 
according to the equations of motion. Solutions to the equations of motion typically approach 
a vacuum configuration uniformly throughout space at late times, and these solutions are 
no exception. There are, however, two qualitatively different possible outcomes of this 
evolution. One possibility is that at some time during the evolution, a, the magnitude of 
the Higgs field, is zero at some point in space and at late times the configuration tends 
toward a vacuum with Nh = 0. As we have seen in Section II, fermion number is not 
violated if the configuration unwinds in this fashion. The second possibility is that at late 
times the configuration dissipates toward a vacuum with Nh unchanged from its initial value 
of ±1. In other words, the Higgs field does not unwind but instead the gauge field "winds 
up". We saw in Section II that fermions are produced in the background of a solution 
which evolves in this way. Without CP violation, for every Nh = +1 configuration which 
relaxes in a baryon producing fashion there is an Nh = — 1 configuration which produces 
anti-baryons. With the inclusion of the CP violating operator O, the hope is that there will 
be some configurations which produce baryons whose CP conjugate configurations relax to 
the Nh = vacuum without violating baryon number. 

We now pause to say a few words about the dynamical context in which we wish to 
use the scenario just described. We imagine that the (thin) bubble wall has just passed, 
leaving in its wake the configuration we are looking at, but that this configuration has 
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not yet had time to relax to equilibrium. Our goal is a qualitative understanding of the 
dynamics of this relaxation. A first order electroweak phase transition can be characterized 
by the change in the gauge invariant quantity (cr 2 ). If we renormalize (a 2 ) such that it is 
equal to v 2 at zero temperature, then for a strongly first order phase transition it is close 
to v 2 just below T c in the low temperature phase, and is much smaller just above T c in the 
high temperature phase. This is a slight motivation for considering initial configurations 
in which a = v throughout space, even though this is in reality not a good description of 
the non-equilibrium configurations left in the wake of the wall and is in fact not maintained 
during the subsequent evolution. There is no justification for choosing either a spherically 
symmetric configuration, or one with = 0, or one which is initially at rest. We shall see 
that this is the Achilles heel of the whole approach. 

It is worth noting that Turok and Zadrozny introduced the analysis we are describing in 
a somewhat different context. They considered the evolution of configurations during the 
passage of a thick wall, rather than after the passage of a thin wall. Instead of using the 
operator ( |1.2| ) to introduce CV violation, they began with the two Higgs doublet model and 
used the operator 

OTiFF (3.1) 

where 6 is a CV odd phase between the two Higgs fields. (Although the operator ( |3.1|) is CV 
even, CV violating dynamics in the Higgs sector picks a sign for the change in (9) during the 
phase transition, thus communicating the CV violation to the gauge sector.) To avoid the 
complications of simulating all the fields in the two Higgs doublet model, they simplified the 
problem as follows. First, they considered circumstances in which (6) was spatially uniform 
but changes linearly in time during the passage of the bubble wall. That is, they replaced 
#(x, t) in ( p.l[ ) by 9(t) = ct with c a constant. Second, they chose to consider the effect of this 
forcing term, now treated as externally imposed, on just the standard model fields. That is, 
they found solutions to the equations of motion obtained from a Lagrangian given by Q2.1p 
with the addition of the term (ctTrFF). Since this Lagrangian depends explicitly on time, 
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there is no conserved energy in their problem. If we wanted to explore the effects of the CV 
violating operator ([L2|) used in this paper as (a 2 ) changes during the passage of the wall, 
we could follow precisely the same strategy. Instead, we wish to treat the non-equilibrium 
conditions after the passage of a thin wall. Therefore, we simply use a Lagrangian given by 
( j2.ip plus ( |1.2| ). This means that energy is conserved (to better than half a percent in our 
numerical simulations) during the evolution of the gauge and Higgs fields. 

We wish to solve the equations of motion obtained from the action ( |2.1|) augmented by 



the addition of the operator ( |1.2|) . We work in the spherical ansatz |30| in which all gauge 



invariant quantities are functions only of r and t, and solve the equations numerically. In 
this paper, we do not describe our numerical methods in any detail; they are standard and 
are similar to those of Ref. [29[|. We also do not present our solutions in full. Rather, we 



describe them only in sufficient detail to demonstrate that it is very difficult to follow this 
approach to completion. Let us first consider initial conditions of the form ( |2.11| ) with r)(r) 
of ( |2.12| ) given by 



r 

1 — tanh ( — 



(3.2) 



where R is a constant parametrizing the size of the configuration. This configuration satisfies 



the boundary condition (|2/7|) and has Higgs winding number Nh = +1, and we use it as 
the initial condition for the equations of motion, setting all time derivatives to zero at 
t = 0. For the moment, we set b = in ( |1.2|) and do not introduce CV violation. In 
agreement with Turok and Zadrozny, we find that there is a critical value of R, which we 
call defined as follows. For all R < R^ the configuration evolves toward a vacuum 
configuration with Nh = 0. No fermions are produced in this background. For all R > R£, 
the configuration evolves toward a vacuum configuration with Nh = +1, and fermions are 
produced. The values of R^ which we obtain in simulations with several different values 
of mn/mw are in quantitative agreement with those obtained in Ref. ||. Repeating this 
exercise beginning with i](r) given by —1 times that in ( |3.2| ), that is beginning with the CV 
conjugate configuration having Nh = —1, we find an analogously defined R~ . As we have 
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not yet introduced any CV violation, we necessarily find R~ = R+. We now repeat the entire 
procedure with b 7^ 0, that is with CV violation present. Although at this point quantitative 
comparison with Ref. || is no longer possible because of the differences described above, we 
do in fact find R~ 7^ R£. Unfortunately, this is not the end of the story. 

We now investigate slightly more general initial conditions. Consider initial configura- 
tions exactly as above except that the time derivative of o is nonzero and is given by 

a{r) = 7^ 2 [1 - tanh (r/i?)] , (3.3) 

with 7 some constant. We now find that for some values of 7, R~ < R£ whereas for other 
values of 7, R~ > Rf. This dooms an analysis in terms of the single parameter R. Clearly, 
a more general framework is needed. 

Consider a family of initial configurations with Nh = +1, much more general than we 
have considered to this point, parametrized by a set of parameters To this point, we 
have introduced two such parameters, R of (|3.2|) and (|3.3j ) and 7 of ( [Ol) . More generally, 
we must consider other profile functions for 77 and & and must allow for initial configurations 
in which a varies with r and and are nonzero. More generally still, one would 
want to go beyond the spherical ansatz. Of course eventually one wants to work towards 
an analysis involving an infinite set of /3's, but it seems reasonable to start with some 
finite set For any fixed b, define a function F + (/3i, /3 2 , ■ ■ •) which has the following 
properties. F + (/3i,/? 2 , . . .) > for all points in /3-space which describe configurations which 
evolve toward the N H = +1 vacuum thereby producing fermions, and F + {(3\ 1 /3 2 , ■ ■ •) < 
for all points in /5-space describing configurations which evolve towards the Nh = 
vacuum. If we only consider initial configurations described by fl2.11|) , (|2.12|) and ( |3.2|) 
which are parametrized by the single parameter R, then we can take F + (R) = R — R+. 
Completely analogously, we define a function F~(/3i, /3 2 , • • •) such that the hypersurface 
F~(Pi,p 2 , . . .) = divides the /3-space of Nh = —1 configurations into those which evolve 
towards the Nh = — 1 and Nh = vacua. 

In the absence of CV violation, /3 2 , . . .) = and F~(f3 1 , /3 2 , . . .) = define the same 

hypersurface. In this case, imagine allowing a CV symmetric ensemble of configurations with 
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(a) (b) 

FIG. 1. Sketch of two qualitatively different possible behaviors of the critical surfaces F + = 
(represented by a solid line) and F~ = (represented by a dashed line) in the space (/3i,^2, ■ ■ •) 
describing initial configurations. 

Nh = +1 and Njj = — 1 to evolve. (By CP symmetric we mean that the probability for 
finding a particular Nh = +1 configuration in the ensemble is equal to that for finding its 
CP conjugate Nh = — 1 configuration.) Since the configurations which anomalously produce 
fermions are exactly balanced by those which anomalously produce antifermions, the net 
fermion number produced after relaxation would be zero. We wish to investigate the behavior 
in the presence of the CP violating term O of ( |1.2| ). The hope is that O will affect the 
dynamics of Nh = +1 configurations and Nh = — 1 configurations in qualitatively different 
ways and that after relaxation to vacuum a net fermion number will result, even though 
the initial ensemble of configurations was CP symmetric. With 6 ^ we have seen that the 
two hypersurfaces F + = and F~ = are indeed distinct. The configurations represented 
by points in /5-space between the two hypersurfaces yield a net baryon asymmetry. There 
are two qualitatively different possibilities, however, which we illustrate schematically in 
Figure 1. In Figure la, the hypersurfaces F + = and F~ = do not cross. The sign 
of b has been chosen such that we produce net baryons in the region between the two 
surfaces. In Figure lb, the hypersurfaces cross and we produce net baryons in regions B 
and D, and net antibaryons in regions A and C. If the hypersurfaces do not cross, as in 
Figure la, then a simple estimate of the fraction of configurations which yield a net baryon 
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asymmetry is possible. This fraction would be proportional to the separation between the 
two hypersurfaces measured in any direction in /3-space with a component perpendicular to 
the hypersurfaces — for example, it would be proportional to — R~) — and it would 
be proportional to b/M 2 , the coefficient of O. Unfortunately, we have seen that when we 
consider the two parameter space of (R, 7) the hypersurfaces F + = and F~ = do in fact 
cross. Thus, in the more general space (fa, (3%, . . .) the picture cannot look like that sketched 
in Figure la and must look like that sketched in Figure lb. 

There are other indications that life is complicated, as in Figure lb. We have found 
that the evolution of configurations near the F + = and F~ = surfaces is in some cases 
extremely complicated. For example, there are instances in which just to one side of F + = 0, 
the configuration changes its winding number from Nh = 1 to Nh = to Nh = 1 to Nh = 
to Nh = 1, having a zero of the Higgs field at four different times, before finally settling 
toward the Nh = 1 vacuum, whereas just on the other side of F + = the configuration 
goes through three zeroes of the Higgs field before settling towards the Nh = vacuum. 
This behavior, also seen by Turok and Zadrozny, suggests that the location of the F + = 
surface in /3-space depends on details of a complicated nonlinear evolution, and so must be 
a highly nontrivial function of many of the parameters specifying the initial configuration. 
It would be very nice to find a simple characterization of the F + = and F~ = surfaces 
in terms of only one or a few parameters, but this seems very unlikely. Hence, although 
the dynamics of the unwinding of topological configurations after the phase transition in 
the presence of the operator O may lead to a baryon asymmetry, we see no way to make a 
simple analytical estimate of this asymmetry. This may well be a valid way of looking at 
the microphysics of electroweak baryogenesis, but it seems that large scale 3 + 1 dimensional 
numerical simulations of the kind recently pioneered by Moore and Turok |31| (but including 
CV violation via ( |1 . 2|) and working in a setting in which the bubble walls are thin and rapidly 
moving) are required in order to estimate the contribution to the BAU. In Section V, we will 
return to a brief discussion of the large scale numerical simulations which seem necessary. 
Before that, in the next section, we explore a different attempt at obtaining a semi- analytical 
estimate of the magnitude of the effect. 
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IV. KICKING CONFIGURATIONS ACROSS THE BARRIER 



Let us now turn to what may seem initially to be a somewhat orthogonal discussion of 
the physics of local electroweak baryogenesis, following that of Dine, Huet, Singleton, and 
Susskind 0. In a sense, this discussion is more general than that of the previous section, 
because it attempts to treat baryon number violating processes of a type more general than 
the unwinding of winding number one configurations. On the other hand, the treatment of 
these more general processes is, of necessity, greatly over-simplified. The complications we 
encountered in the previous section are real, but they seem not to appear in the treatment 
of this section, and we are able to obtain an estimate for the baryon asymmetry generated. 
A priori this could either mean that the method extracts the essence of the matter or that 
it sweeps important physics under the rug. The truth is somewhere in between these two 
extremes. We give a discussion of the physics under the rug, and argue that the estimate 
we obtain should be seen as an upper bound on the baryon asymmetry generated by local 
electroweak baryogenesis in the thin wall scenario. Although we follow Ref. |7| to some 
extent, our discussion does not exactly parallel theirs and we shall note the points where we 
differ as we come to them. 

In the high temperature phase, baryon number violating processes are not exponentially 
suppressed. The barrier crossing configurations typically [H|P2"| | have sizes given by the 
magnetic correlation length 

^(awT)- 1 . (4.1) 

We will think of dividing space up into cells of this size, and looking at configurations cell 
by cell. The energy in gauge field oscillations with wavelength £ is of order T, but the total 
energy in a cell is much larger, as it is presumably of order T 4 £ 3 . Indeed, this energy is much 
larger than the sphaleron energy, which is i£ S ph ~ v/g. Most of the energy is in oscillations of 
the gauge and Higgs fields on length scales shorter than £. These configurations are crossing 
the barrier between vacua via regions of the barrier far above the lowest point on the barrier, 
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that is far above the sphaleron, and they look nothing like the sphaleron. It was initially 
thought JO] | that in each cell of volume £ 3 , the sphaleron barrier was crossed once per time 



£, leading to a baryon number violation rate per unit volume of the form (|1 . 1|) with k ~ 1, in 



agreement with numerical simulations |33|]. A recent analysis []32| in fact suggests that the 
time it takes for a configuration in a given cell to cross the barrier is of order ti/otw leading 
to k ~ aw- We now consider what happens when the bubble wall hits the configurations 
just described. 

Consider the configuration in one cell. It traverses a path through configuration space, 
which we parametrize by r. Dine et al. consider the special case in which this path is the path 
in configuration space which an instanton follows as a function of Euclidean time r, but this 
is not essential, and it is clear that they were thinking of more general circumstances also. 
The configurations discussed in Section III can be seen as special cases of those described 
here. The energy of the configuration has a maximum at some r (at which the configuration 
crosses the barrier) which we define to be r = 0. Following Dine et ai, we now write down a 
Lagrangian which is intended to describe the dynamics of r as a function of time for r near 
r = 0: 

C( T , f) = 2 + -^-r 2 + C -^—n 2 r . (4.2) 
2£ 2£ 3 £ M 2 K J 

In this expression, Cj, C2, and C3 are dimensionless constants, different for each of the infinitely 
many possible barrier crossing trajectories. The factors of £ have been put in by dimensional 
analysis treating r as a quantity of dimension —1. (We shall see that rescaling r by a 
dimensionful constant does not change the final result.) This Lagrangian should be seen as 
the first few terms in an expansion in powers of r and f. Because we have assumed that 
t = is a maximum of the energy as a function of r, no odd powers of r can appear. In the 
absence of CP violation, there can be no odd powers of f, since they make the dynamics 
for crossing the barrier from left to right different than from right to left. The TtFF in 
the operator O includes a term which is proportional to the time derivative of the Chern- 
Simons number, and this means that O must contribute a term in L which is linear in f. It 
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is obviously quite an over-simplification to treat barrier crossing as problem with one degree 
of freedom. As we saw in Section III, the complete dynamics can be very complicated, even 
for relatively simple initial conditions. We will return to this point, but for now we forge 
ahead with (gg). 

The momentum conjugate to r is given by 

P = J T + JJp a (4 ' 3) 
and the Hamiltonian density is therefore 

«^HM-# r2 - (4 - 4) 

Before considering the thin wall case of interest in this paper, it is worth pausing to consider 
the thick wall limit in which (a 2 ) is changing slowly and other quantities evolve adiabatically 
in this slowly changing background. A reasonable assumption is that the variables (r, p) are 
Boltzmann distributed with respect to the Hamiltonian (|4.4|) at each instant, treating a 2 as 
approximately constant. This implies that the distribution of p is centered at 

However, from (|4.3|) we see that this means that the velocity f is Boltzmann distributed with 
center f = 0. Thus, the presence of the CP violating operator ( |1.2|) does not bias the velocity 
of trajectories in configuration space in the thick wall limit. This conclusion disagrees with 
that of Ref. 0. There is nevertheless an effect. Integrating the third term in (f4.2|) by 
parts, one obtains a term linear in r proportional to the time derivative of a 2 . This changes 
the shape of the potential energy surface in configuration space during the passage of the 
wall, and yields an asymmetry. In this limit, in which the wall is thick and departure from 
equilibrium is small, the problem is much more easily treated in the language of spontaneous 
baryogenesis 0,03 — the operator O acts like a chemical potential for baryon number. 

We now turn to the thin wall case. Immediately after the wall strikes, the fields are not 
yet in equilibrium. As we saw in the previous section, life is complicated. The idea of this 
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section, however, is to use an impulse approximation to estimate the kick which f receives as 
the wall passes, and from this to estimate the baryon asymmetry that results. The equation 
of motion for r obtained from ( |4. 2| ) is 

c 2 c 3 b d 2 . . 

During the passage of a thin wall, the first term on the right hand side can be neglected 
relative to the second. Making this impulse approximation, we find that the passage of the 
wall kicks f by an amount 

A, = , (4.7) 

where Act 2 is the amount by which a 2 changes at the phase transition. The kick Af has a 
definite sign. Thus, in the thin- wall limit, the distribution of the velocities in configuration 
space of barrier crossing trajectories is biased and a baryon asymmetry, whose magnitude 
we now discuss, results. 

If Af is large compared to f , the velocity the configuration would have had as it crossed 
t = in the absence of the action of the wall, then Af will kick the configuration over the 
barrier in the direction it favors, and will produce, say, baryons rather than anti-baryons. If 
Af is small compared to to, it will have no qualitative effect. The fraction of the distribution 
of configurations with f < Af is proportional to Af . (Note that had we taken r to have 
dimension other than —1 in ( |4.2j ), various powers of £ would have run through the calculation 
until this point when they would have cancelled in computing the fraction of the distribution 
of f affected by the wall.) Note that in this calculation it was not necessary for r to be 
precisely at r = when the wall hits. It was only necessary for r to be close enough to 
t = that the Lagrangian ( |4.2j ) is a good approximation. It is difficult to quantify what 
fraction of configurations satisfy this criterion of being "close enough to r = 0" , so we simply 
parametrize our ignorance by calling this fraction /, for fudge factor. It is worth noting that 
/ does not depend on the time it takes configurations to traverse the barrier. As discussed 
at the beginning of the section, this is now thought to be of order £/acw- Making this time, 
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say, longer just makes the time during which the configuration is "close enough to r = 0" 
longer and need not affect /. Tabling further discussion of / momentarily, we estimate the 
net number density of baryons produced as 



where we have absorbed the constants q into /. 

At the time of the electroweak phase transition, the entropy density of the universe is 
s ~ 45T 3 , and so we obtainQ 



to increasing a 2 up to that value at which baryon number violating processes become expo- 
nentially suppressed in thermal equilibrium. Whereas in the thick wall case, baryon number 
violating processes stop when a 2 reaches this value, this is not the case in the thin wall 
scenario. In this setting, thermal equilibrium is not maintained even approximately, and 
we see from the above discussion that what matters is the net change in a 2 as the wall 
passes. (If we assume that a 2 changes arbitrarily rapidly as we did by using the impulse 
approximation above, then the final answer can only depend on the total jump in a 2 and 
cannot depend on the value of a 2 at which equilibrium baryon number violation ceases.) 
Once one picks an extension of the standard model which makes the transition strongly first 
order, one can compute Act 2 . Here, we will simply take Act 2 = v 2 /2, which is approximately 
what is obtained in the minimal standard model with a 35 GeV Higgs mass |21| . Putting it 
all together, we find 



n B ~Atf£, 



-3 



(4.8) 




(4.9) 



The size of the effect clearly depends on Aa 2 . It has been suggested that Act 2 corresponds 




(4.10) 



1 This result agrees with that of Ref. 0, although our discussion and theirs are somewhat different. 
In comparing to Ref. |7| note that they take b ~ aw- 
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If we take, for example, b ~ aw and M ~ 1 TeV, the bound ( |1.5| ) can be satisfied and (f4. 10|) 
suggests that a cosmologically relevant BAU may be generated. We see that if CV violation 
is introduced via the operator O with a coefficient b/M 2 satisfying ( |1.5| ), and if the bubble 
walls are thin, then the contribution to the baryon asymmetry of the universe from local 
electroweak baryogenesis can be at an interesting level so long as the fudge factor / is not 
smaller than about a tenth. 

We have reached the estimate ( |4.10|) by arguing that tib must be proportional to 
(b/M 2 )Aa 2 , arguing by dimensional analysis that it must be proportional to £~ 3 , argu- 
ing that the time it takes configurations to cross the barrier in equilibrium in the high 
temperature phase is not relevant, and lumping our remaining ignorance into /. We now 
turn to a discussion of what goes into /, and how the treatments of Sections III and IV 
are related. There are many contributions to /, since the treatment leading to the estimate 
Q4.10Q is greatly over-simplified. First, there are the constants q, which of course differ for 
the different configurations in different cells of volume £ 3 , and must somehow be averaged 
over. Second, using the impulse approximation is not really justified. In reality, the wall 
does not have zero thickness. More important, even if the wall is thin, the time during which 
it can affect a configuration of size £ is at least £. Third, as we have already mentioned, the 
treatment in terms of the Lagrangian ( j4.2|) only has a chance of capturing the physics near 
t = 0, and it is not at all clear what fraction of configurations satisfy this. Configurations 
which happen to be farther away from the crest of the ridge between vacua when the wall 
hits do receive a kick from the wall. However, even if this kick is large, it may not be in 
a suitable direction in configuration space to be effective. Configurations far from r = 
can contribute to rig, but their contribution is hard to compute, because there is no way to 
reduce the problem to one of one degree of freedom far from r = 0. Fourth, even near the 
crest of the ridge for a given trajectory the problem does not really reduce to one degree 
of freedom. For the configurations of interest, a is a function of space and time and the 
operator O and the bubble wall conspire to affect its dynamics. We have attempted to 
describe the effect by treating a as constant in space and time on either side of the wall and 
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only changing at the wall. This is a caricature at best. Fifth, by now the problem should 
be sounding more like that of Section III, and we must face up to the specific difficulties 
discussed there. After the wall has passed, the fields are not yet in thermal equilibrium and 
their dynamics is complicated. This may in fact yield a further contribution to rig. It may 
also, however, negate some of the contribution estimated in ( 4.10|) because some configu- 



rations kicked across the barrier in one direction by the passage of the wall may at a later 
time wander back across the barrier whence they came. As we discovered in Section III, an 
estimate of the magnitude of these sorts of effects is difficult even for the restricted class of 
configurations we considered there. To sum up, / is almost certainly less than 1. Hence, it 
would be best to use ( 4.10|) as an upper bound on ub/s, rather than as an estimate. 



V. CONCLUDING REMARKS 

If the electroweak phase transition is strongly first order, it is possible that the observed 
baryon asymmetry of the universe may have been generated by non-equilibrium processes 
occurring as bubble walls sweep through the plasma. In general, there will be contributions 
to the asymmetry both from local baryogenesis and from nonlocal baryogenesis. In the last 
few years, much effort has been devoted to computing the contribution from nonlocal effects 
PJ9|-|14||. Recently, the work of Moore and Turok |JT| has made it clear that the computa- 



tional resources now available make a large scale numerical treatment of local baryogenesis 
a possibility. Our goal in this paper was to reconsider two possible routes to a semi-analytic 
estimate of the magnitude of the effect. If the bubble walls are thick, conditions remain close 
to thermal equilibrium during the passage of the wall, and the nonequilibrium physics can 
be captured by assigning nonzero chemical potentials to various quantum numbers includ- 
ing baryon number. Analytic estimates for the BAU produced in this setting exist in the 
literature J7|J| and the need for a numerical treatment is not pressing. If the bubble walls 
are thin, however, or if (as is no doubt the case) they are comparable in thickness to other 
length scales in the problem, the situation is unsettled. We have reanalyzed the problem of 
local baryogenesis in the thin wall limit using two different semi-analytical approaches 
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The method of Dine et al. does yield an estimate ( 4.10|) , but the uncertainties parametrized 



by /, particularly those highlighted by the difficulties which prevented us from obtaining an 
estimate via the method of Turok and Zadrozny, mean that this should be viewed at best 
as an upper bound. A large scale numerical treatment seems necessary. 

Moore and Turok |JT| have recently taken a big step in this direction. They have per- 



formed 3 + 1 dimensional classical simulations in which a bubble wall moves through a box 
converting the high temperature phase to the low temperature phase. To date, they have 
focused more on computing quantities like the wall thickness, the wall velocity, the surface 
tension, and the drag on the wall and have only begun their treatment of local electroweak 
baryogenesis. To this point, they have introduced CV violation only by "mocking up" the 
effects of ( |1.2j ) by first computing the average wall profile (cr)(z) for an ensemble of walls, 
and then doing a simulation in which one measures the distance of a given point to the 
nearest bubble wall and adds a chemical potential for Chern-Simons number at that point 
proportional to the spatial derivative of the average wall profile at that distance. This chem- 
ical potential is only nonzero on the wall, as it would be if it were proportional to j|cr 2 for a 
moving wall. Nevertheless, by imposing the chemical potential as an external driving force 
instead of simply introducing ( |1.2| ) in the Lagrangian and letting the dynamics do their thing 
self-consistently, one risks missing a lot of the difficulties (and potential effects) we have dis- 
cussed in Section III and at the end of Section IV. The simulations of Moore and Turok 
suggest that a large scale numerical assault on the problem of local electroweak baryogenesis 
is now possible; the difficulties we have discussed which prevent us from obtaining a reliable 
semi-analytic estimate of the magnitude of the effect show that it is necessary. 



Although we have taken considerable care to describe the failings of ( 4.10|) , even if we use 
it only as an upper bound it is interesting in the following sense. Let us assume that whatever 
CV violation is introduced in order to make baryogenesis possible can be parametrized 
by the CV violating operator O of ( |1.2| ). Combined with the experimental bound ( |1.5| ) 
on the coefficient of O, the result ( |4.1U| ) shows that if the experimental sensitivity to the 
electric dipole moment of the electron or the neutron can be improved by about an order 
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of magnitude, and if these experiments continue to yield results consistent with zero, then 
the baryon asymmetry of the universe produced by local electroweak baryogenesis is smaller 
than that observed, even if future numerical simulations were to demonstrate that / is as 
large as 1. 
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